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Abstract. We report on two different manipulation procedures of a tunable rf 
SQUID. First, we operate this system as a flux qubit, where the coherent evolution 
between the two flux states is induced by a rapid change of the energy potential, turning 
it from a double well into a single well. The measured coherent Larmor-like oscillation 
of the retrapping probability in one of the wells has a frequency ranging from 6 to 20 
GHz, with a theoretically expected upper limit of 40 GHz. Furthermore, here we also 
report a manipulation of the same device as a phase qubit. In the phase regime, the 
manipulation of the energy states is realized by applying a resonant microwave drive. 
In spite of the conceptual difference between these two manipulation procedures, the 
measured decay times of Larmor oscillation and microwave-driven Rabi oscillation are 
rather similar. Due to the higher frequency of the Larmor oscillations, the microwave- 
free qubit manipulation allows for much faster coherent operations. 
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1. Introduction 

Superconducting qubits are promising systems for the realization of quantum 
computation. Their coherent evolution, entanglement, storage and transfer of quantum 
information as well as quantum non-demolition readout are just few examples of what 
has been already achieved [U EJ [3l Hj. One of the major challenges superconducting 
qubits are facing now is to increase their coherence time in order to reach an useful 
number of quantum operations [5]. However, there is an alternative possibility of 
achieving the same goal: one can try to increase the qubit operation speed to make 
qubit gates shorter, thus leading to less restrictive requirements for the coherence time. 
Nowadays, qubit operation frequencies exceeding 500 MHz are difficult to reach with 
phase qubits manipulated by microwave signals, since the limited anharmonicity of the 
qubit energy potential gives rise to a leakage into higher excited states. 

In this paper, we report experiments with a tunable rf SQUID and demonstrate 
two different manipulation procedures. The system under test, a double SQUID, can be 
operated in two different regimes: (i) as a flux qubit manipulated via fast (dc) pulses 
of magnetic flux, and (ii) as a phase qubit in which the quantum state evolution is 
controlled by microwave pulses of chosen amplitude and phase. The former manipulation 
approach without using microwaves allows to reach very high oscillation frequencies, 
while operating the same circuit as a phase qubit offers the possibility to verify the 
obtained results using a well defined and studied manipulation technique. We will see 
in the following that the coherence times measured using the two different manipulation 
schemes are rather similar, suggesting that decoherence acts in a similar way in both 
cases. This conclusion draws attention to common sources of decoherence (presumably, 
dielectric loss due to two-levels fluctuators in the material used for chip fabrication 
[6l [7]), and emphasizes the relevance of increasing the number of qubit rotations within 
the coherence time by increasing the oscillation frequency of the system. 

The paper is organized as follows. A detailed overview of the system under test, 
composed of the double SQUID and the readout dc SQUID, is given in section [2j 
Section [3] describes the manipulation of the double SQUID by deforming its energy 
potential using fast pulses of magnetic flux. The measured coherent oscillations are 
presented here together with their theoretical interpretation. Section H] reports the 
manipulation of the same device as a phase qubit. Rabi oscillations measured via 
microwave pulses of variable length are reported. The dependence of the Rabi oscillation 
frequency on microwave power is analyzed by taking into account the population of 
higher excited states. 

2. System and experimental setup 

The circuit that we studied consists of a superconducting Nb loop of inductance L = 85 
pH interrupted by a small dc SQUID of inductance I = 6 pH. The name for this kind 
of device varies between different authors from 'double SQUID' [8J to 'combined rf-dc 
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SQUID' [9 J or 'modified rf SQUID' [10J. The two Josephson junctions embedded into the 
dc SQUID are nominally identical except for unavoidable asymmetries originating in the 
fabrication process. Each junction has a critical current of 8 /zA and a self capacitance of 
0.4 pF. The system is manipulated via two magnetic bias fluxes $ x and $ c applied to the 
large and small loops, respectively. The detailed characterization of the device taking 
into account non-identical Josephson junctions and non- negligible inductance I has been 
reported elsewhere [llj. Both loops are designed fully gradiometrically with the intent 
to decrease both the noise induced by external uniform magnetic fields as well as cross 
talk between the two bias fluxes. A photograph of the double SQUID together with 
its schematic representation is presented in figure (Ha). The double SQUID is defined 
by the pale-white area delimited by white solid lines at the center of the picture. The 
two large holes of 100 x 100 pm 2 define the gradiometric main loop. In the pale-white 
area highlighted by the white dashed ellipse one can identify two much smaller holes 
(10 x 10 /xm 2 ), defining the gradiometric inner dc SQUID. The two Josephson junctions, 
of dimensions 3x3 /im 2 , are visible at the center. The two coils on the left and right 
side are used to control the bias flux <& x \ the mutual inductance between them and the 
double SQUID is 2.6 pH. The coil inducing the bias flux $ c , visible on the lower-central 
part of the picture, is wrapped around one of the two small loops, to which it has a 
mutual inductance of 6.3 pH. The placement of coils and holes reduces the cross talk 
between <3> x and $ c lines to less than 1%, as verified experimentally. Each of the two 
coils inside the large loops forms a part of a superconducting transformer connecting the 
system to an unshunted readout dc SQUID (one for each side). The circuit was made 
by Hypres [12] using standard Nb/AlO^/Nb technology, with a critical current density 
of 100 A/cm 2 and SiC>2 as dielectric material for junction isolation. 

The inductance of the small loop I is chosen to be much smaller than that of the 
main loop L, so that the two-dimensional energy potential defining the system can 
be approximated by a one-dimensional function of the parameter S, corresponding to 
the phase difference across the inner dc SQUID. Moreover, the inductance of the main 
loop L and the critical current of each junction Jo, are chosen such that the parameter 
Pl — 2vrL ■ 2J /$o ranges between 1 and 5tt/2, resulting in a double well potential. The 
two bias fluxes & x and $ c are used to manipulate the energy potential profile. Changes 
on $> x modify the symmetry (figure DJb)) of the potential, while changes on <3> c tune the 
height of the barrier between the two local minima (figure QJc)). 

Measurements were performed at the base temperature of a dilution fridge stabilized 
at 15 niK. The experimental setup is reported in figure [2j The two currents generating 
the fluxes $> x and $ c were supplied via coaxial cables with 10 dB attenuators anchored at 
the lK-pot stage of the refrigerator. The microwave signal, applied to the double SQUID 
via the coil <& x , was supplied through an extra coaxial cable with 10 dB attenuation at 
the lK-pot stage and 10 dB at base temperature. For biasing the readout dc SQUID we 
used superconducting wires, metal powder filters [13] at base temperature and current 
dividers at the lK-pot stage, while the lines to measure the voltage response were 
equipped with low pass filters instead of current dividers at the lK-pot stage. 
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Figure 1. (a) Photo and schematic representation of the studied double SQUID. The 
entire device occupies a space of 430 x 230 fim 2 . (b)-(c) Manipulation of the energy 
potential profile via the two bias fluxes <5> x and <£><,. Variation of the symmetry is 
achieved by changes of the flux $> x , while the barrier between the two local minima is 
tuned by the flux $ c . 
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Figure 2. Representation of the experimental setup. The pale-blue area indicates the 
device on chip, while the symbols Pw-fil, divider, RC filter and 10 dB indicate powder 
filters, current dividers, RC filters and attenuation values on coax lines. 
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Figure 3. Variation of the potential shape in the microwave- free manipulation of the 
double SQUID. 

3. Microwave-free manipulation 

The qubit manipulating scheme reported in this section is based on changing the double 
well potential to a single well shape and back, caused by fast dc pulses on the bias flux 
$ c . Hereby, the computational states are mapped to the flux states of the double well 
potential (i.e., left and right wells), while the coherent evolution occurs between the two 
lowest energy eigenstates of the single well potential. 

The manipulation procedure is described in detail in Ref. [H]. It consists in the 
four main steps depicted in figure [3](c) and outlined as follows. 

® The state of the qubit is initialized. A pulse on the flux Q x tilts the potential 
in order to remove one of the two minima. This strongly asymmetric potential is 
maintained for a time required for the complete relaxation in the energy minimum. 

© After the initial qubit state is prepared, the potential is turned into the symmetric 
double well. The qubit state remains frozen due to the large barrier height between 
the wells, which prevents any tunneling between them. In this case, only pure 
dephasing is possible as only one of the two states is populated. Which of the two 
wells is populated depends on the sign of the flux pulse in the previous step. 

(3) Here we perform the quantum manipulation: By applying a short flux pulse $ c 
the barrier is completely removed and the potential is turned from the symmetric 
double well into single well. This deep single well can be approximated by a 
parabola of a harmonic oscillator, whose characteristic frequency depends on the 
pulse height. When the pulse is over, the potential returns to the initial deep 
double-well configuration. 

(§) The final circuit state is read out by applying a current ramp through the dc SQUID 
and recording its switching current to the non-zero voltage state. 

The main point of above procedure is the non-adiabatic transition between the 
lowest two energy levels, which occurs during the potential transformation. In figure H] 
the energy values, expressed as a frequency, of the first six eigenvalues for the double- 
SQUID are reported versus the control flux $ c . On that graph the energy of the ground 
state |0) is defined to be the zero reference level. Symmetric double well situation and 
single well situation are depicted, respectively, on the left and right side of the graph. 

For a perfectly symmetric double well, the lowest two energy levels are degenerate, 
with an energy splitting A between them depending exponentially on the inverse of 
the barrier height. For the more realistic case of non perfect symmetry, we express 
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Figure 4. Energy values of the first six eigenvalues for the double-SQUID. The ground 
state |0) is defined to be the zero reference level. 



the asymmetry e as the energy difference between the two minima of the potential. For 
large potential barrier, the level splitting hfi, = \/e 2 + A 2 is dominated by e and remains 
essentially constant during the initial drop of the barrier. 

When the barrier is almost completely removed (for e ~ A, which is indicated 
in figure H] as "portal" region [15J), the level splitting starts to increase and saturates 
at the single- well oscillation frequency. This process corresponds to "half a Landau- 
Zener process, starting from the degeneracy point and arriving at the large energy gap 
condition (see figure H]). In this process, the system is prepared and remains initially 
in a 50% superposition of the two lowest energy eigenstates. The initial phase of this 
superposition depends on the chosen initial state of the deep double-well potential. For 
example, an initial left \L) (right \R)) flux state is transformed in a 50% superposition 
with positive (negative) sign. Although unwanted non-adiabatic transitions to upper 
levels (from the third onward) are possible, they can be suppressed owing to the gap 
existing between the first two levels and upper levels in the portal region. This requires 
an appropriate choice of the pulse rise/fall-times so that the transition is performed 
non-adiabatically for the computational state but adiabatically for upper levels. 

After crossing the portal region, the system is maintained for the pulse duration in 
the deep single-well potential. Here, the lowest two levels are equally populated and the 
relative phase between them evolves with a rate given by the level splitting. This level 
splitting is tunable by the flux pulse amplitude, and can be described by the formula 



where /3($ c ) = Pl cos(7r$ c /$ ) is the modification of the parameter (3l via the flux $ c . 
During this time the system is weakly responsive to fluctuations on the flux $ c and so it 
is naturally protected against noisq§|. The final phase between the states is determined 

§ as an example, a fluctuation of l/i$o on the flux <I> C is responsible of a change in the oscillation 
frequency of only 125 kHz, corresponding to a percentage variation of 6.6 x 10~ 4 at 19 GHz. 




(1) 
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Pulse width (ns) 

Figure 5. Probability of measuring the state \L) as a function of the pulse duration. 
The coherent oscillation shown here has a frequency of 14 GHz and a coherence time 
of approximately 1.2 ns. 

by the duration of the flux pulse, i.e. by the time elapsed in the single well condition. 
At the end of the pulse the flux returns to the initial condition and the system goes back 
to the two well state. The portal region is crossed again, and the inverse of the previous 
process occurs: the relative phase between the two energy states is transformed to the 
amplitude of the left/right states. After the pulse, the system's flux state becomes once 
again frozen, and it is read out by an unshunted dc SQUID magnetometer which is 
weakly inductively coupled to the qubit loop. 

The described procedure was repeated many times (from 100 to 10000) in order 
to determine the probability of one of the projected states for a chosen combination of 
parameters, such as pulse height and duration. This is then repeated for different 
pulse durations At in order to record the coherent oscillations of the qubit state 
(figure [5]). From curves collected for different pulse amplitudes we extracted the 
oscillation frequency and the decay time. In figure [6] we show the measured oscillation 
frequency (full dots) versus the pulse amplitude. These data are in very good agreement 
with the theoretical expectation given by ([T|) for the single- well oscillator frequency (solid 
line). The measured decay time is on the order of a nanosecond, independent of the 
frequency. 

The described manipulation scheme only allows for ^-rotations around the qubit 's 
Bloch sphere. The full qubit control requires also phase control achieved by R z - 
rotations. These, in principle, can be realized by slightly unbalancing the potential 
through a dc pulse on the bias flux $> x . Quantum computation protocols require also 
the interaction between different qubits which can be envisaged by using an inductive 
switchable coupling that allows for o z o z coupling between them |16j . 

It is worth noting that the first manipulation of a superconducting qubit with fast 
dc pulses was achieved by Nakamura and collaborators on a Cooper-pair box in 1999 
[17] . In that article, coherent oscillations of frequencies up to 19 GHz are reported. In 
spite of the similar manipulation technique the coherent oscillations reported by us are 
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Figure 6. Measured oscillation frequencies versus amplitude of the short flux pulse 
(full dots). The solid curve is a numerical simulation using the measured parameters 
of the circuit. 

due to a different physical process. In the charge qubit, the manipulation was done 
at the degeneracy point where ground and first excite state, defined by n and n + 1 
cooper pairs on the island, are no eigenstates of the energy. The population of the 
initially prepared ground state is thus oscillating with the first excited state due to the 
off-diagonal terms in the Hamiltonian. The coherent oscillations reported in [T7j are 
a direct measurement of the population evolution. In our system the free evolution is 
performed in a single well potential where there are no degeneracy points. The phase 
evolution of each populated energy level is the only physical process taking place. The 
phase difference between ground and first excite state is not detectable in the energy 
base, but it can be easily measured by projecting the state of the system in the flux 
base defined by the double well situation. The reported coherent oscillations are thus a 
direct consequence of the not directly measurable phase evolution in the energy base. 

4. Microwave-induced manipulation 

In addition the above microwave-free manipulation, we have also operated the system 
as a phase qubit using microwave driving. The operating procedure is the same as 
that reported in Ref. [18]. The energy potential profile is strongly tilted via the flux 
bias <& x making one of the two local minima shallow enough to contain only a small 
number of energy levels. The two computational states of the qubit are defined to 
be ground |0) and first excited states |1) in the shallow well. The first excited state 
is populated by resonant absorption of photons from the microwave field, while the 
complete manipulation on the Bloch sphere can be performed via microwave pulses of 
defined duration and phase and dc flux pulses. The main difference between our device 
and the conventional rf SQUID phase qubits [19] is the possibility to tune in situ the 
Josephson energy of the device via the flux coil $ c . This additional tuning parameter 
allows to modify the anharmonicity of the energy potential in a slightly different way 
to what is done by the flux Q x . It is thus possible to minimize the leakage to higher 
excited states in a strongly deformed potential (3-4 energy levels in the shallow well) 
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Figure 7. Rabi oscillation of the double SQUID manipulated as a phase qubit by 
applying microwave pulses at 19 GHz. The oscillation frequency changes from 540 
MHz to 1.2 GHz by increasing the power of the microwave signal by 10 dB. 



without changing the number of levels inside the well nor the escape probability of each 
state. The readout procedure is performed in two steps: 

(i) An adiabatic, but fast, dc pulse is sent to the bias coil $ x , the potential is thus 
deformed and consequently the barrier separating the two wells is reduced. The 
amplitude of the pulse is calibrated such that the transition from the shallow to the 
deep well is triggered for the ground state |0) with a probability of approximately 
10%. That ensure a complete escape of the system from the first excited state |1), 
leading to a theoretical visibility close to 90%. 

(ii) The flux of the double SQUID is measured with the dc SQUID coupled inductively 
to it, using the procedure already described in section [3j 

The sequence of preparation, manipulation and readout is repeated between 100 
to 10000 times, depending on the desired statistical error. An example of three Rabi 
oscillations obtained for different powers of the microwave driving field at 19 GHz is 
reported in figure [7J The increase of the oscillation frequency with microwave power is 
clearly seen. The topmost oscillation has a frequency of /R a bi — 540 MHz. Oscillations 
at lower frequencies are difficult to measure since the decay time is very short, between 
1.5 and 2.0 ns. On the lower graph, the oscillation frequency is approximately equal to 
/Rabi = 1-2 GHz. Further increasing the power leads, as we will see, to an unwanted 
population of higher excited states. 

The correlation of the Rabi oscillation frequency versus amplitude of the microwave 
driving field is reported in figure El The deviation from the linearity visible at higher 
powers is a clear indication of populating higher excited states [20] . We note also that the 
qubit has have been driven slightly off-resonance since the frequency of Rabi oscillations 
does not reach zero for zero power. 

The experimental data in figure M were fitted with two curves: 
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Figure 8. Measured Rabi oscillation frequency versus the normalized amplitude of the 
microwave signal (solid dots). The dashed line is a linear fit taking into account slightly 
off-resonance microwave field, while the fit represented by the solid line considers a 
population of higher excited states. 



• The dashed line represents a linear power-frequency relation with the addition of a 
term taking into account the microwave pumping out of resonance. The equation 
has the form 

/Rabi = Vt/Ol - /Ml) 2 + (kAf 

where A is the normalized amplitude of the microwave driving signal in mW 1 / 2 
and k is a constant. The microwave frequency /mw — 19 GHz is off set from 
the resonance frequency /oi corresponding to the transition between the lowest 
two levels. The points used for the fit refer to a normalized microwave amplitude 
ranging between 0.1 mW 1 ^ 2 and 0.4 mW 1 ^ 2 . From the fit we obtained a value of 
k = (2919 ± 105) MHz/mW 1 / 2 and f 01 - f MW = (349 ± 45) MHz. 

• The solid fitting curve takes into account both a non-zero population of higher 
excited states due to large driving field and the manipulation red-detuned from the 
resonance. The dependence of the Rabi frequency on the power of the microwave 
signal is described e.g. in Ref. [21J (equation 22). The function that we used for 
the solid curve fit has the form 

/Rabi = V%01-/MW) 2 + (£^(1-/M 2 )) 2 . 

As before, (/ i — /mw) is the frequency out of resonance and k the linear power- 
frequency correlation, while the new coefficient (3 depends on the details of the 
system. The result of the fit, this time made through all measured points, gives 
the values of k = (2977 ± 104) MHz/mW 1/2 , f 01 - f MW = (350 ± 42) MHz and 
(3 = 0.411 ± 0.060 mW" 1 . Within the noted errors, both fits yield the same values 
for the common parameters. 

We measured the energy relaxation time T\ of the double SQUID when it was 
operated as a phase qubit. The occupation probability of the first excited state was 
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Figure 9. Measurement of the relaxation time T\ for the double SQUID operated as 
a phase qubit. 

measured after a variable time between a resonant microwave 7r-pulse and the readout 
pulse. The measured probability decay was fitted exponentially yielding T\ = 1.37 ns 
(figured]). 

5. Conclusion 

We presented measurements on a double SQUID manipulated both as a double- well flux 
qubit and as a phase qubit. The device manipulation as a flux qubit was performed by 
modification of its energy potential profile via fast dc pulses, while the manipulation as 
a phase qubit required the use of microwave signals to induce Rabi oscillations between 
the ground and first excited state of a shallow well in a strongly deformed potential. 
The measured coherence time of the Larmor oscillations obtained in the flux regime is 
about 1 ns, of the same order as the relaxation time T\ in the phase regime. Since best 
available phase qubits already display relaxation times of several 100 ns, obtained by 
using appropriate materials in the fabrication processes [6], we suppose that also the 
coherence time of the Larmor oscillations obtained without microwaves could strongly 
benefit from the same treatment. Such a possible improvement of the coherence by 
two orders of magnitude, together with the much higher oscillation frequency in the 
microwave-free Larmor mode, should in principle allow to reach the ultimate goal of 
10 4 single-qubit gate operations within the coherence time which is needed for the 
implementation of quantum algorithms. 
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